Resonating-valence-bond structure of Gutzwiller-projected superconducting wave 
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Gutzwiller-projected (GP) wave functions have been widely used for describing spin-liquid physics 
in frustrated magnets and in high-temperature superconductors. Such wave functions are known to 
represent states of the resonating-valence-bond (RVB) type. In the present work I discuss the RVB 
structure of a GP singlet superconducting state with nodes in the spectrum. The resulting state for 
the undoped spin system may be described in terms of the "path integral" over loop coverings of the 
lattice, thus extending the known construction for RVB states. The problem of topological order in 
GP states may be reformulated in terms of the statistical behavior of loops. The simple example of 
the projected d-wave state on the square lattice demonstrates that the statistical behavior of loops 
is renormalized in a nontrivial manner by the projection. 
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Gutzwiller-projected (GP) wave functions have been 
often applied to describe the unconventional properties 
of high-temperature superconductors and the spin-liquid 
phase of frustrated spin- 1/2 systems 0, H H 3- Since 
the early days of their use for strongly correlated sys- 
tems, GP wave functions have been known not only to 
provide a good variational ansatz, but also to represent 
the resonating-valence-bond (RVB) physics of the ground 
state. Because of their RVB structure, the GP wave func- 
tions merit deep investigation as independent objects, ir- 
respective of the underlying physical Hamiltonian d Q . 

It was shown in the original anaysis of GP supercon- 
ducting wave functions Q that such wave functions cor- 
respond to the RVB states with singlet amplitudes ay 
given by the Fourier transform of 



a(k) 



Uk_ 

Vk 



(1) 



where Uk and Vk are the coherence factors of the BCS 
wave function before projection. For a fully gapped su- 
perconductor, this procedure produces amplitudes ay 
decaying exponentially with the site separation \i — j\, 
and the interpretation as a RVB state is straightfor- 
ward. However, in the case of a BCS superconductor 
with nodes, including the commonly used d x 2_ y 2 state, 
this simple derivation produces values of a(k) singular at 
the nodal points, and the resulting singlet amplitudes ay 
are long-ranged. Therefore, in the case of a supercon- 
ductor with nodes, the conventional derivation cannot 
interpret the GP wave function as a RVB state. 

In this note I resolve the problem of interpreting the 
physical content of the GP wave function by using an al- 
ternative representation in terms of the "loop-soup" path 
integral. In this construction, well-known for RVB states, 
the correlation functions are expressed as sums over all 
close-packed coverings of the lattice with closed loops 
01 In the conventional RVB loop construction, 
the statistical weights are determined by the products 
of singlet amplitudes along the loops. In my generalized 
derivation, the product of singlet amplitudes is replaced 
by the trace of the product of the BCS Green's functions 



along the loop. Thus the role of the RVB singlet ampli- 
tudes is played by the BCS Green's functions (2x2 ma- 
trices). In the case of a superconductor with nodes, the 
BCS Green's functions are only marginally local: they 
decay algebraically with distance. It is therefore clear 
that caution is required when applying the spin-liquid 
RVB scenario to high-temperature superconductors [lOj. 

The paper is continued with a brief discussion of the 
relation of the loop-soup construction to the topological 
order predicted for RVB states H El El E3 The 
presence of topological order corresponds to disordered 
short-ranged loops. 

I then employ the example of the undoped GP <i-wave 
state on the square lattice to demonstrate the possibility 
of strong renormalization of the loop statistics by pro- 
jection. The loop-length distribution may be character- 
ized by spin-spin correlations which decay algebraically 
with distance. The power of the algebraic decay is non- 
universal and depends on the variational parameter of 
the wave function. 

The derivation begins by denning the GP wave func- 
tion. I consider a system consisting of spins 1/2 occupy- 
ing a finite lattice with an even number of sites N. The 
GP wave function is constructed by using an auxiliary 
Hamiltonian of the BCS form, 



(2) 



(in this general form the hopping and pairing amplitudes 
between any pairs of sites are arbitrary). The Hamilto- 
nian is spin-rotation invariant (it involves only singlet 
pairing), and its ground state ^bcs is a spin singlet. 
The GP wave function 'J'gp for the undoped system is 
constructed by projecting onto states with exactly one 
fcrmion per lattice site (Gutzwiller projection). The re- 
sulting wave function may be written as a wave function 
of spin variables ^cp({ai}), where the spins Oi must con- 
tain N/2 up and N/2 down spins. 
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The Hamiltonian (0) contains more information than 
required for the construction of the wave fucntion. First, 
the Gutzwiller projection selects states invariant with re- 
spect to SU(2) gauge rotations in the particle-hole space 
0,0. Different Hamiltonians related by such SU(2) 
gauge symmetries therefore produce the same projected 
wave function ^gp- Second, the Hamiltonian con- 
tains information not only about its eigcnfunctions, but 
also about the spectrum, which is not used in the con- 
struction of the wavefunction. To eliminate the latter re- 
dundancy, we shall use the equal-time Green's function 
(the "projector operator") instead of the Hamiltonian. 

To construct the equal-time Green's functions, we hrst 
define the 2iV-dimensional vector space of fermionic op- 
erators 



E 



[u(i)Vl T + v(i)ipix 



(3) 



(the sum is taken over the N lattice sites). This space 
has the Hermitian form defined by the anti-commutator 
{7i:72j-, and the Hamiltonian acts in this space as 
a Hermitian operator. The eigenvalue equations are the 
Bogoliubov-de Gennes equations 



[H n ] = Ej. 



(4) 



The spectrum of the eigenvalues consists of pairs of op- 
posite energies ±E: this is the consequence of the spin- 
rotational invariance of the Hamiltonian J5J . For each so- 
lution (u, v) with energy E, the pair (v*, —u*) gives a so- 
lution with energy —E. Thus there are N positive-energy 
solutions and N negative-energy solutions to Eq. @. 

Next we define the projector onto the negative-energy 
states, 



G= J2 fa> fa 



(5) 



£ fc <o 



where \jk) are normalized eigenstates [solutions of 
The operator G constructed in this way is a 2N x 2N 
Hermitian matrix which may be considered as a single- 
particle projector onto the negative-energy states or, 
equivalently, as the matrix of equal-time Green's func- 
tions. It contains no information about the eigenval- 
ues of QJ, but only about the eigenfunctions. Thus the 
projector G contains less information than the original 
Hamiltonian, but remains sufficient for constructing the 
multi-particle ground-state wave function. 

In real space, the operator G may be represented by 
the set of 2 x 2 matrices Gij for each pair of lattice sites 
(i,j). The hermiticity implies 



(6) 



The SU(2) gauge transformation acts on G by conju- 
gation 



» w}G l3 W h 



(7) 



where Wi are the SU(2) matrices of the gauge transfor- 
mation. 

From the completeness of the basis of all states 7& 
and from the E — > —E symmetry (i.e. from the spin- 
rotational invariance), we can derive the following iden- 
tity for G^: 



G i: j + J'G*jJ = Sijl, 



J = 





-1 



(8) 



Mathematically, this identity expresses the fact that G — 
2 belongs to the Lie algebra sp(iV) (see, for example, 
Ref. Ilfih . One may verify that Eq. JHJ is invariant with 
respect to SU(2) gauge rotations. 

From the SU(2) invariance of Gutzwiller projection, 
it follows that the wave function depends only on the 
SU(2)-gauge- invariant properties of the Green's functions 
Gij . For future application, we define the cyclic trace of 
the Green's functions 



^i 1 ...ii — l ^T(Gi 1 i 2 Gi 2 i 3 



■ ^Mi) 



(9) 



which depends on the oriented closed loop . . . ij. This 
type of cyclic trace appears below in the loop path inte- 
gral. 

To proceed with the loop construction, we first express 
the wave function in terms of the negative-energy states 
7fc and further develop the loop path integral for ground- 
state expectation values. 

The BCS ground state of the Hamiltonian may be 
expressed as 



*bcs = n u -■ i) 



(10) 



E k <0 



where |J. ... |) is the state with all lattice sites occupied 
by down spins, || ... j) = Yl % iph |0), and the operators 
7fc are the negative-energy eigenvectors of Eq. Q . 

After projecting onto singly-occupied states, this wave 
function may be shown to yield for the GP wave function 



*cp({c r J) = Det fc ,i [u k (i)\v k (i)] 



(11) 



where [uk(i)\vk{i)) is the N x N matrix whose first N/2 
columns are composed of the coefficients Uk(i) while the 
last N/2 columns involve those of ujt(i). The index k 
referes to all negative-energy states and the index i to all 
spin-up sites in the given spin configuration {en}. 

Alternatively, the same wave function may be written 
in the same form 1)11(1 , but with the index i labeling spin- 
down sites. It follows from the rotation invariance that 
these two determinants must give the same wave func- 
tion, to within an overall phase factor. 

The experssion for the "partition function" {\& | \T/} is 
obtained by taking for the bra-vector the determinant 
(|ll|l over the spin-up sites and for the ket- vector the same 
determinant over the spin-down sites. The product of the 
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FIG. 1: (a) An example of covering the lattice by loops in the 
"loop-soup" construction, (b) To each link of the loop, there 
corresponds a 2 x 2 matrix Green's function Gy. Each (ori- 
ented) loop is assigned the amplitude —2Ti 1 ...i l . The origin of 
the negative sign in the loop amplitude lies in the fermionic 
statistics in the GP construction. The vertices are labeled 
with up or down spins which alternate along the loop. The 
factor of two in the loop amplitude is due to the two possible 
spin labelings of the loop. 

two determinants is then rewritten as the determinant of 
the product of the two matrices (taking the sum over the 
index k) to yield 

|*gp({^})| 2 = const DetijiGij], (12) 

where i labels spin-up sites and j labels spin-down sites. 
The matrix is the N x N matrix composed of 2 x 2 
blocks Gij. The constant prefactor does not depend on 
the spin configuration {<7i} and will be omitted in fur- 
ther calculations. The expression (|12fl is explicitly SU(2) 
gauge- invariant . 

The determinant l|12|) may be further expanded as a 
product of non-intersecting loops and summed over all 
possible spin configurations {<7i} to produce the partition 
function (\f f |\f f ). Brief algebraic manipulations using the 
symmetry (JSJ give the loop path integral 

<*!*>= En(- 2T <i-o> ( i3 ) 

{C„} n 

where the sum is taken over all coverings by non- 
intersecting loops of even length, and Ty...^ are the loop 
traces defined in Eq. © , see Fig. ^ The sum (|13fl [and 
all loop sums below] is to be understood as that over 
oriented loops. From the symmetry (jHJ it follows that 
the traces T^...^, for even- length loops are purely real 
and do not depend on the orientation of the loop. How- 
ever, it is important to preserve loop orientations in the 
definition of the sum (|13|) . because they affect the mul- 
tiplicities of length- two loops: every length- two loop has 
only one orientation and appears in the sum (|13fl once; 
all longer loops admit two orientations and appear twice. 

This result has the form of a path integral in which dif- 
ferent correlation functions may be computed. For exam- 
ple, the spin-spin correlation function (S z (i)S z (j)) may 
be written as 

{C*„} n 

(14) 




FIG. 2: Defining topological sectors on a cylinder. The "ref- 
erence line" (dashed) connects the two edges of the cylinder, 
(a) The "plus-minus" sectors in the loop path integral are 
generated by changing the sign of Green's functions Gy in- 
tersecting the reference line, (b) The "even-odd" sectors in 
the RVB states are selected as those with an even (odd, re- 
spectively) number of singlets intersecting the reference line. 

where the sum is now taken only over loop coverings with 
the sites i and j belonging to the same loop, and (— l) Pij 
takes values ±1 depending on whether the number of 
loop links between i and j is even or odd. 

The path integral (|13|l has the same form as the loop 
construction for RVB states. Most generally, the loop 
path integral may be written as 

(*i*)= En^«)' ( i5 ) 

{C„} n 

where A{C n ) is an amplitude depending on the geometry 
of the loop C n . 

This general formulation includes many different RVB- 
type wave functions in a variety of systems. The GP wave 
functions in Anderson's derivation [lj correspond to the 
loop amplitude 

A(Cn) — —2a,i 1 i 2 ai 2 i 3 . . . (Xijjj , (16) 

where the singlet amplitudes defined from ljT|l . 

Note the importance of the negative sign in the loop 
amplitudes (|13fl and l|16fl . This sign arises from the 
fermionic statistics involved in the GP construction. Al- 
ternatively, for a spin system it is possible to construct a 
"bosonic" product of singlets which leads to a loop am- 
plitude similar to Eq. 1161) , but without the negative sign 

Remarkably, the Rokhsar-Kivelson ground state of 
dimer models ^lj may also be described with the same 
formalism l|15|) : it is sufficient to set A(C n ) equal to 1 
for length-two loops on allowed dimer positions and to 
otherwise. 

The concept of the "loop soup" aids in visualizing the 
conditions required for the topological order proposed for 
RVB states [ll[ll[3. As explained in Refs.HH the 
different topological sectors on a multiply connected do- 
main (e.g. on a cylinder or a torus) may be accessed by 
imposing periodic or antiperiodic boundary conditions on 
the fermions in the Hamiltonian J5| along a topologically 
nontrivial contour. As a result, the periodic/antiperiodic 
boundary conditions produce two different GP wave func- 
tions ^+ and >F_. The topological order implies two dual 



4 




FIG. 3: Main panel: staggered spin correlations in the un- 
doped projected d-wave state for different values of the varia- 
tional parameter <f> — 4 arctan A. The Monte Carlo simulation 
was performed on the 24x24 lattice and involved averaging 
over 10 4 samples. The error bars are smaller than the sym- 
bol sizes. Inset: The exponent a as defined in Eq. <|23J as 
a function of the variational parameter (f>. Squares: values 
of a obtained from the spatial decay of spin correlations in 
the 24 x 24 system. Stars: values of a obtained from the 
size-dependence of the squared integrated staggered magneti- 
zation M(L) defined in Eq. 11241 . see Fig. 0] below. The error 
bars are of the order of the symbol sizes. The arrow marks 
the position of the optimal value of the variational parameter 
(4> ~ 0.34-7T) minimizing the variational energy of the Heisen- 
ber g H amiltonian within the given class of wave functions 



conditions 

(V+\X\V + ) - -> 0, (17) 

(*+|X|*_)^0, (18) 

where the correlation functions must tend to zero suffi- 
ciently rapidly (e.g., exponentially) with increasing sys- 
tem size for any local operator X. 

It will become clear below that the first of these con- 
ditions (|17fl may be viewed as an effective absence of in- 
finite loops in the path integral l|15|) , and may be related 
loosely to the short range of spin correlations. The sec- 
ond condition (|18fl corresponds to the absence of valence- 
bond crystallization in the RVB construction, and, in the 
general case (|15f) . may in all probability be formulated as 
the absence of loop crystallization. Note, however, that 
the following discussion of those correspondences is not 
fully rigorous and some statements are indicated as con- 
jectures. 

I begin by considering the condition Ijl7|l . Under the 
assumption that that the Green's functions Gij decay 
rapidly with distance, loops of large size may appear only 




L 

FIG. 4: The dependence of the fluctuating integrated stag- 
gered magnetization M(L) defined in Eq. 124H on the linear 
system size L in the log-log scale. The values of a deduced 
from these data are plotted by star symbols in the inset of 

Fig. m 

with many links of short length. In this case, the change 
of boundary conditions in Eq. I|13|l may be represented by 
changing the sign of the Green functions Gij intersecting 
the reference line (Fig. EJ) . Consequently, the partition 
functions in the ± sectors may be written as 

(*±i*±)= ^iirnn^-O' (19) 

where W is the total winding number of the loops. This 
definition may be extended directly to the general for- 
malism (|15|) and to any expectation value of a local ob- 
servable X, 

^ ± \X\^ ± )=J2(^) W X({C n })l[A(C n ). (20) 

{c„} 

Therefore the difference between the two topological sec- 
tors Ijl7(l contains only configurations with odd winding 
numbers. The absence of large loops then serves as a 
sufficient condition for the criterion (|17J) of topological 
order [I believe that in most situations this is also a nec- 
essary condition]. On the other hand, the absence of 
large loops also guarantees that spin correlations (|14fl are 
short-ranged, but is not a necessary condition. Therefore 
one may expect that in many cases the short range of spin 
correlations is related to the condition l|17fl . but in this 
contribution I provide no rigorous derivation of such a 
relation. 

We consider next the second condition l|18|) . For the 
RVB wave functions represented as linear combinations 
of products of singlets, this condition corresponds to the 
absence of singlet ordering (e.g. the absence of a valence- 
bond crystal (l7J). This can be shown by employing the 
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even-odd basis denned conveniently in the case of RVB 
or dimer states on a topologically non-trivial domain [l3l 
Hij . Using the expression (|2*U|) for the definition of \P±, 
the "plus-minus" basis is related to the "even-odd" basis 
by 



(21) 



where ^ and ^f e are the states with even and odd num- 
bers of intersections of the singlets (dimers) with the ref- 
erence line (Fig. |2J). The condition l|18|l may then be 
reformulated as 



(* e |X|* e ) - (* |X|* ) -»0 



(22) 



which, in turn, can be interpreted as the absence of sin- 
glet (dimer) crystallization [l8| . One may expect that for 
the more general loop construction l|15|l this condition 
also corresponds to the absence of loop crystallization. 
However I do not have a rigorous argument supporting 
this conjecture. 

As discussed above, the first criterion of the topologi- 
cal order (|17J) involves the question of the effective loop 
size in the path-integral ensemble l|15|l . This effective 
loop size cannot be inferred easily from the singlet am- 
plitudes ciij (in the case of the RVB construction) or from 
the Green's functions Gij . Numerical studies of spin cor- 
relations in "bosonic" RVB states Q indicate that the 
effective loop size may be renormalized nontrivially by 
the constraint of fully-packed loops. Similar renormal- 
izations of loop statistics may also be observed in the 
fermionic GP construction. I consider for illustration 
one of the simplest examples of a GP wave function, the 
undoped projected d-wave state on the square lattice. 
This state has been studied previously as an approximate 
ground state of the Heisenberg antiferromagnet 0, H(| . 
It does not possess long-ranged antiferromagnetic order, 
but only a power-law decay of antiferromagnetic correla- 
tions ja, l21| . The only dimensionless parameter in the 
wave function is the pairing strength A — A/i, and 
the wave function is invariant under A — > A -1 (see, for 
example, Ref. fl9|). Equivalently, the same wave func- 
tion may be defined as the projected normal staggered- 
flux state with a flux per plaquette of <fi = 4 arctan A 
[l"! Efi . l22j | . Numerical Monte Carlo simulations indicate 
that for all values of the parameter <p, the antiferromag- 
netic correlations obey the power law, 



(-iy-i(S z (i)S*(j))K\i-j\ 



(23) 



see Fig. [3] One observes that the spin correlations decay 
very slowly (with a < 2), which implies that the effective 
loop size in the loop path integral l|15l) is of the order 
of the system size. The condition ()17|l of the topological 
order is therefore expected to fail. Interestingly, in Ref. 
this state was indeed classified as one without topo- 
logical order, but on the basis of the failure of the other 



condition (|18fl [the condition (|17fl was not considered for 
this state]. Note that in the current example (undoped 
projected d-wave state) the spin correlations l|23")) probe 
precisely the probability of the two sites i and j belong- 
ing to the same loop: the Green's functions Gy connect 
only sites of opposite antiferromagnetic sublattices, and 
therefore the term (— l) Pij in 1|14|) is equal to (— l) 4- - 7 
independently of the loop configuration. 

The actual value of the exponent a in the power law 
(|23|1 remains a subject of controversy. One possible sce- 
nario is the universal value of a depending only on the 
symmetries of the wave function (for example, the au- 
thors of Ref. propose the value a = 3/2). However, 
my numerical results suggest a different possibility: a 
parameter-dependent non- universal exponent a(0). The 
values a(<fi) extracted from the spatial decay of spin cor- 
relations are shown in the inset of Fig. EI To verify those 
results, I have also computed the square of the integrated 
staggered magnetization in the L x L system, 



(24) 



The dependence of M{L) on the linear system size L is 
plotted in Fig. 0] for three values of <f>. A simple scaling 
argument predicts M(L) oc L 2 ~ a . The values of a(4>) 
deduced from M(L) are plotted by star symbols in the 
inset of Fig. While those results are consistent with a 
</>-dependent value of a, they do not constitute a definite 
proof: a more detailed analysis of finite-size effects or an 
analytic argument is needed to settle the issue. 

To summarize, the GP superconducting wave functions 
admit a loop description generalizing the conventional 
loop construction for RVB states. The role of the sin- 
glet amplitudes is played by the equal-time BCS Green's 
functions. The properties of topological order and frac- 
tionalization are related to the correlations of the loops, 
but these cannot be inferred simply from the unprojected 
Green's functions involved. Nevertheless, the loop formu- 
lation may be helpful in visualizing the RVB properties 
of the GP wave functions, and possibly for developing 
appropriate analytic approximations. 

The author thanks G. Jackeli and B. Normand for help- 
ful comments on the manuscript. 
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